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Abstract 

Research on stabilization of coupled hyperbolic PDEs has been dominated by the focus on pairs of 
counter-convecting (“heterodirectional”) transport PDEs with distributed local coupling and with controls 
at one or both boundaries. A recent extension allows stabilization using only one control for a system 
containing an arbitrary number of coupled transport PDEs that convect at different speeds against the 
direction of the PDE whose boundary is actuated. In this paper we present a solution to the fully general 
case, in which the number of PDEs in either direction is arbitrary, and where actuation is applied on 
only one boundary (to all the PDEs that convect downstream from that boundary). To solve this general 
problem, we solve, as a special case, the problem of control of coupled “homodirectional” hyperbolic 
linear PDEs, where multiple transport PDEs convect in the same direction with arbitrary local coupling. 
Our approach is based on PDE backstepping and yields solutions to stabilization, by both full-state and 
observer-based output feedback, trajectory planning, and trajectory tracking problems. 
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1. Introduction 

a) Background: Coupled first-order linear hyperbolic systems, typically formulated on a 
1-D spatial domain normalized to the interval (0,1), are common in modeling of traffic flow 0, 
heat exchangers fl2^ . open channel flow DU, S or multiphase flow DTOl . [fT4l . ifTSll . 

Research on stabilization of such PDEs has been dominated by the focus on pairs of counter- 
convecting transport PDEs with distributed local coupling. In DU a first solution allowing actua¬ 
tion on only one boundary and permitting coupling coefficients of arbitrary size was presented. 
A recent extension IfT^ by three of the authors of the present paper allows stabilization using 
only one control for a system containing an arbitrary number of coupled transport PDEs that 
convect at different speeds against the direction of the PDE whose boundary is actuated. 

In this paper we present a solution to the fully general case of coupled hyperbolic PDEs. We 
divide such PDE systems into two categories: 

. homodirectional systems of m transport PDEs, for which all the m transport velocities have 
the same signs, i.e., all of the PDEs convect in the same direction. Because of the finite 
length of the spatial domain, these are inherently stable but the coupling between states can 
cause undesirable transient behaviors and the trajectory planning problem is non-trivial. 

. heterodirectional systems of n -l- m transport PDEs, for which there exist at least two 
transport velocities with opposite signs, i.e., where m PDEs convect in one direction and n 
PDEs convect in the opposite direction. The coupling between states traveling in opposite 
directions may cause instability. 

In this paper we present control designs for the fully general case of coupled heterodirectional 
hyperbolic PDEs, allowing the numbers m and n of PDEs in either direction to be arbitrary, and 
with actuation applied on only one boundary (to all the m PDEs that convect downstream from 
that boundary). To solve this general problem, we solve, as a special case, the heretofore unsolved 
problem of control of coupled homodirectional hyperbolic linear PDEs, where multiple transport 
PDEs convect in the same direction, have possibly distinct speeds, and arbitrary local coupling. 

Our approach is based on PDE backstopping and yields solutions to stabilization, by both full- 
state and observer-based output feedback, trajectory planning, and trajectory tracking problems. 

b) Literature: Controllability of hyperbolic systems has first been investigated using explicit 
computation of the solution along the characteristic curves in the framework of norm DT^ . 
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uni, Bill. Later, the so-ealled Control Lyapunov Funetions methods emerged, enabling the design 
of dissipative boundary eonditions for nonlinear hyperbolie systems in the eontext of both 
norm and norm IH, BH, [O. Further, using Lyapunov funetions method, suffieient boundary 
eonditions for the exponential stability of linear [fT3l or nonlinear [fTTl . [fT8l hyperbolie systems 
of balanee laws have been derived. All of these results impose restrietions on the magnitude of 
the eoupling eoeffieients, whieh are responsible for potential instabilities. 

In [|71, a full-state feedbaek eontrol law, with aetuation only on one end of the domain, whieh 
aehieves exponential stability of elosed-loop 2-state heterodireetional linear and quasilinear 
hyperbolie systems is derived using a baekstepping method. With a similar baekstepping trans¬ 
formation, an output-feedbaek eontroller is designed in [fT^ for heterodireetional systems with 
m = 1 (eontrolled) negative veloeity and n (arbitrary) positive ones. These results hold regardless 
of the (bounded) magnitude of the eoupling eoeffieients. Unfortunately, the method presented in 
0, m ean not be extended to the ease m > 1. 

c) Contribution: The first step towards this paper’s general solution for m > 1 was presented 
(but not published as a paper) in Ii24l for m = 2 and n = 0. In eonferenee paper [[20ll . an extension 
to m = 2 and n = 1 is aehieved. 

The eontribution of this artiele is two-fold. For (n -l- m)-state heterodireetional systems, we 
derive a stabilizing boundary feedbaek law that ensures finite-time eonvergenee of all the states 
to zero. For homodireetional systems (for whieh stability is not an issue), we design a boundary 
eontrol law ensuring traeking of a given referenee trajeetory at the uneontrolled boundary. 

Both designs rely on the baekstepping approaeh. A partieular ehoiee of the target system, 
featuring a easeade strueture similar to [|71 Seetion 3.5], enables the use of a elassieal Volterra 
integral transformation. Well-posedness of the system of kernel equations, whieh is the main 
teehnieal ehallenge of this paper, is proved by a method of sueeessive approximations using a 
novel reeursive bound. 

In the ease of heterodireetional systems, the approaeh yields a full-state feedbaek law that 
would neeessitate full distributed measurements to be implemented, whieh is not realistie in 
praetiee. For this reason, we derive an observer relying on measurements of the states at a single 
boundary (the anti-eontrolled one). Along with the full-state feedbaek law, this yields an output 
feedbaek eontroller amenable to implementation. 
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d) Organization: In Section]^ we introduce the model equations. In Sectionwe present 
the stabilization result for heterodirectional systems: the target system is presented in Sec¬ 
tion |III-A| while the backstepping transformation is derived in Section |III-B[ The design is 


summarized in Section III-C In Section IV we present the boundary observer design. In Sec¬ 


tion |V] we present the motion planning result for homodirectional systems. Section VI contains 
the main technical difficulty of the paper, i.e. the proof of well-posedness of the backstepping 
transformation. We conclude in Section |V^ by discussing open problems. 

II. System description 

We consider the following general linear hyperbolic system 

Ut{t, x) + A^Uxit, x) = x) + x) 

Vt{t, x) - A~Vx(t, x) = ir^u{t, x) + 'L~~v{t, x) 
with the following boundary conditions 

u{t, 0) = Qov(t, 0), v(t, 1) = Riu{t, 1) -I- U(t) 

where 

u = [ui ■ ■■ Un] , V = Vi ■ • ■ r 
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Remark 1: We consider here constant coupling coefficients and transport velocities for the 
sake of readability. The method straightforwardly extends to spatially varying coefficients, with 
more involved technical developments. 

Besides, we also make the following assumption without loss of generality 

Vj=l,...,m o-~~ = 0, (11) 

i.e. there are no (internal) diagonal coupling terms for v-system. Such coupling terms can be 
removed using a change of coordinates as presented in, e.g., [|7]| and Il20ll . This yields spatially- 
varying coupling terms, which is not an issue in the light of Remark 

Remark 2: If two or more states have the same transport speeds (i.e. /r, = nj for some i j) 
we refer to those states as isotachic. This case was intentionally avoided in Q. To deal with 
isotachic states, we consider the change of coordinates v(r, x) = A(x)v{t, x). The matrix A(x) is 
a block-diagonal matrix, with An = 1 if /r, jxj for j 4^ i. If there is a set of n, isotachic states 
(i.e. there is i such that /r, = /i, for y = / -l-1,..., / -l- n, - I, then there is in A{x) a corresponding 
block B{x) of dimension n, x n, in A{x). Each of these B{x) is computed independently for each 
isotachic set of states. If we call 2,^0 the matrix of coupling coefficients among these isotachic 
states (i.e. with coefficients cr” for j,k = i,i + I,... ,i + rii - 1), then B{x) is computed from the 
initial value problem B{x) = I!niBixyLiso, B{0) = . It is easy to see that this transformation 

is invertible, since one can define a matrix C(x) from C(x) = I jHi'LisoBix), C(0) = /„,xn, . One has 
that C{x) is the inverse of B{x) as B(0)C(0) = /„,x«, and ^B(x)C(x) = 0. Applying this invertible 
transformation eliminates the coupling coefficients between isotachic states, but results in some 
spatially-varying coupling terms, which is not an issue as explained in Remark 

III. Stabilization of heterodirectional systems 

In this section, we derive a stabilizing feedback law for the general (n + m)-state system. 
Notice that this is interesting only in the case n 4 0, since instability arises from coupling 
between states traveling in opposite directions. Following the backstepping approach, we seek 
to map system Q-Q to a target system with desirable stability properties using an invertible 
Volterra transformation. 
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A. Target system 

1) Target system design: We map system Q-Q to the following target system 

at{t, x) + IsTax{t, x) = TT^a{t, x) + 'Z^~j3(t, x) 

+ f C^(x,Oa(Od^+ f C-(x,mOd^ 

Jo Jo 

M,^)-A-j8x(t,x) = G(x)m 


with the following boundary eonditions 

a ( t , 0 ) = Qoj 8 ( t , 0 ), 

where and C~ are L“ matrix funetions on the domain 


j 8 ( t , 1) = 0 


r = {o<^<x<i}, 

while G 6 L“(0,1) is a lower triangular matrix with the following structures 


( 12 ) 

(13) 

(14) 

(15) 


G(x) 


0 

g 2 , l ( x ) 




0 




(16) 


The coefficients of C and G will be determined in section TTT-B 


2 ) Stability of the target system: The following lemma asseses the finite-time stability of the 
target system. 

Lemma 3.1: Consider system ([T^,([T3]) with boundary conditions ( pA] ). Its zero equilibrium is 
reached in finite time t = tf, where 


tf 


1 1 

:=l + yl 

bj 


(17) 


Proof: Noting ( |13| )-([T4l) and we find that the jS-system is in fact a cascade system, 
which allows us to explicitely solve it by recursion as follows. The explicit solution of (Si is 
given by 


l-x 


J3i(0,x + /iiit) ift< 
Mbx)=> "" 

10 if t 

Ml 


(18) 
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Notice in particular that ySi is identically zero for t From the time t > on, we have 

that yS 2 (t, x) satisfies the following equation 

/32((t, X) - ^l2p2x{t, x) = 0. (19) 


Similarly, by expressing the solution along the characteristic lines, one obtains that 


J32(t,x) = 0 

Thus, by mathematical induction, one can easily get that /3j(j = 1, • • • ,m) vanishes after 


( 20 ) 


= E-- 


This yields that 


m 

J3it,x) = 0, t>V—. 

7=1 


( 21 ) 


( 22 ) 


When t > Yi the or-system becomes 

at(t, x) + A^axit, x) = S'^'^Q:(t, x) + f C^ix, ^)a(^)d^ 

Jo 

with the boundary conditions 


(23) 


Qf(t,0) = 0. (24) 

Since there are no zero transport velocities for the or-system (see @), we may change the status 


of t and jc, and Equations (23) can be rewritten as 


ax(t, x) + (A^)-^a,(t, x) = x) + I (A^)-^C^(x, ^)Qr(^) 


. +\-lv++ 


f 


L +\“l /^+/ 


(25) 


with the initial condition ( |24l ). Then by the uniqueness of the system (|24l),([25l), and noting the 
order of the transport speeds of the or-system (see (|^), this yields that or identically vanishes 
for 


1 1 
»>-+y- 


(26) 


This concludes the proof. 
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B. Backstepping transformation 

To map system ([T])-@ to the target system (fT^-(fT^, we eonsider the following baekstepping 
(Volterra) transformation 

a{t,x) =u{t,x) (27) 

fit,x)=v(t,x)- f [K(x,mO + L(-^,mf)]df (28) 

Jo 

where the kernels to be determined K and L are defined on the triangular domain 7~. Deriv¬ 
ing ( [28] ) with respect to space and time, plugging into the target system equations and noticing 
that j3(t, 0) = v(t, 0) yields the following system of kernel equations 


0 =K(x, x)A^ + A-K(x, x) + 2“^ 

(29) 

0 =A~L(x, x) - L(x, x)A~ + 

(30) 

0 =K(x, 0)A-^Qo + G(x) - L{x, 0)A“ 

(31) 

0=A-K,(x,f)-K^(x,f)A^ 



(32) 

Q=A-L,(x,^) + L^(x,f)A- 


-L(x,^)2---iS:(x,^)2^- 

(33) 

and yields the following equations for C~(x,f) and C"'(x,f) 


C-(x,0 = L(x,f) + ^ C-(x, s)L(s,f)df 

(34) 

C\x, ^) = K(x, f) + J^ C-(x, s)K(s, 

(35) 


Remark 3: For each x e [0,1], Equation ( |34| ) is a Volterra equation of the second kind on [0, x] 
with C~(x,-) as the unknown. Besides, Equation ( j^ explicitly gives C'^(x,^) as a function 
of C~(x,f) and K(x,f). Therefore, provided the kernels K and L are well-defined and bounded, 
so are and C~. 


Developing equations (|29|)-(|33]) leads to the following set of kernel PDEs 
for I < i < m, 1 < j < n 

n m 

Pid^Kifx, - Ajd^Kijix, = crtjKtkix, f) + J] cr^p]Lip(x, f) 

k=\ p=\ 


(36) 
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for 1 < / < m, 1 < j < m 


liid^Lij{x, ^)+ixjd^Ljjix, ^ o-pj Lip(x, ^ Kik(x, 

p=i 


( 37 ) 


k=\ 


along with the following set of boundary eonditions 


o". ^ 

VI < / < m, 1 < j < n, Kij(x, x) = -^ = ki 


Hi + dy 


‘J 


VI < i,j <m,i + j, 
'i\ < i < i < m. 


Lij(x, x) =- - — = h 


‘J 


Hi - Hj 

n 

HjLij(,X,0) — ^ ^kKjjci^X, O^Cjjc j 


(38) 

(39) 

(40) 


k=l 


To ensure well-posedness of the kernel equations, we add the following artifieial boundary 
eonditions for Lij{i > j) 

L,y(l,^) = lij, for 1 < 7 < / < m (41) 

While the gij, for 1 < j < i < n, are given by 

n 

gij(x) = HjLijix, 0) - ^ ApqpjKipix, 0) (42) 

p=\ 

provided the K and L kernels are properly defined by iUHlI}. whieh we prove in the next 
seetion. 

Remark 4: The ehoiee of imposing ( |4T] ) as the boundary eondition for L,y(l < j < i < m), on 
the boundary x = 1 is arbitrary and was designed to ensure eontinuity of some of the kernels. This 
degree of freedom in the eontrol design had never appeared in previous baekstepping designs 
for hyperbolie system 0, IfT^ . The impaet of the boundary values of L,y, 1 < j < i < m on the 
transient behavior of the elosed-loop system remains an open question, out of the seope of this 
artiele. 

Remark 5: If there are isotaehie states, and the transformation explained in Remark is 
applied, then the kernels for i, j eorresponding to isotaehie states (/i, = Hj) have all boundary 
eonditions of the type (40) instead of (39)—whieh would beeome singular—or ( [4T] ). The results 
that follow do not ehange, but we have omitted the ease for the sake of brevity. 

The well-posedness of the target system equations is assessed in the following Theorem. 
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Theorem 3.2: Consider system ([3^-(|4T]). There exists a unique solution K and L in U°{T). 
Moreover, all the boundary traees for the ^-kernel and L-kemel are funetions of L“(0,1). 

The proof of this Theorem is the main teehnieal diffieulty of the paper and is presented in 
Seetion 


C. Control law and main stabilization result 

We are now ready to state the main stabilization result as follows. 

Theorem 3.3: Consider system ([T])-Q with boundary eonditions Q and the following feed- 
baek eontrol law 


U(t) = -R\u{t, 1) + 


f'm. 

Jo 




(43) 


For any initial condition (uo,Vo) g (T“(0, i))(«+'”)x(«+'”)^ the zero equilibrium is reached in finite 
time t = tf, where tf is given by ( fTT] ). 

Proof: First, notice that evaluating transformation ( [^ at jc = 1 yields ( |43] ). Besides, 
rewriting transformation ([28l) as follows 


a{t, x) 


u{t, x) 

- r' 

V_ 


< 

V_ 

Jo _ 


0 0 
K{x,f) L{x,f) 


u{t, 


d^. 


(44) 


one notices that it is a classical Volterra equation of the second kind. One can check from, 
e-g-, IfT^ that there exists a unique matrix function R. 6 (/^“( 7 “))(”+'”)x(«+m) that 


/ \ 
u{t, x) 


( \ 
a(t, x) 

r 

s_ 

— 

R(x,f) 

, v{t,x) ^ 


, fiit, x) ^ 

Jo 

1 Ab^) ) 


d^. 


(45) 


Applying Lemma |3T] implies that {a, 13) go to zero in finite time t = tp , therefore, by ( |45] ), (m, v) 
also converge to zero in finite time. ■ 


IV. Uncollocated observer design and output feedback controller 

In this section, we derive an observer that relies on the measurement of the v states at the left 
boundary, i.e. 


yit) = v(t, 0) 


(46) 


Then, using the estimates from the observer along with the control law ( [43] ), we derive an output 
feedback controller. 
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A. Observer design 

The observer equations read as follows 

Ut{t, x) + x) x) + x) - P^{x){Kt, 0) - v(t, 0)) (47) 

Vt{t, x) - hTVxit, x) =lr""u{t, x) + lr~v{t, x) - P~{x){v{t, 0) - v(t, 0)) (48) 

with the following boundary eonditions 

u{t, 0) = eov(b 0), v(t, 1) = Riu{t, 1) + u{t) (49) 

where P^{-) and P~{-) have yet to be designed. This yields the following error system 

Ut{t, x) + A^Ujc(t, x) =L'^*u{t, x) + £'^“v(t, x) - P^{x)v{t, 0) (50) 

Vf(t, x) - A~Vx(t, x) ='L~*u{t, x) + £~~v(t, x) - P~(x)v(t, 0) (51) 

with the following boundary eonditions 

M(t, 0) = 0, v(t, 1) = 1) (52) 


Remark 6: One should notice that the output is directly injected at the left boundary, which 
means potential sensor noise is only filtered throughout the spatial domain. Combining the 
approach of IfT^ and the cascade structure of (fT^-(fT4l), we now derive a target system and 
backstepping transformation to design observer gains P""{-) and /’“(•) that yield finite-time 
stability of the error system ([50l)-(|5^. 

B. Target system and backstepping tranformation 

We map system ([50l)-(|52l) to the following target system 

at(t, x) + A^&xit, x) ='L'^'^a(t, x) + f D^{x,^)a{^)d^ (53) 

Jo 

Pt{t, x) - A~^x(t, x) ='L~^a(t, x)+ f D~(x, ^)a{^)d^ (54) 

Jo 

with the following boundary conditions 

a{t,0) = 0, yS(t,l)=i?iS(t,l)- r H(,m^)d^ (55) 

Jo 
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where and D are L“ matrix funetions on the domain T and H 6 L“(0,1) is an upper 
triangular matrix with the following strueture 


H{x) = 


0 h,2W 


0 


0 


(56) 


all of whieh have yet to be determined. 

Proposition 4.1: The solutions of system 
preeisely, one has 


-([T?]) eonverge to zero in finite time. More 


4t > tp. O' = = 0 


(57) 


where tp is defined by ( fTT] ). 


Proof The system eonsists in a easeade of the or-system (that has zero input at the left boundary) 
into the yS-system. Further, the yS is a easeade of its slow states into its fast states. The rigorous 


proof follows the same steps that the proof of Lemma 3.1 and is therefore omitted here. 


To map system (|50l)-(|52|) to the target system (|53])-(|55]), we eonsider the following baekstep- 
ping (Volterra) transformation 


u{t, x) = a{t, x) + f M(x, 

Jo 

-*x 


(58) 

Kt,x)=[5{t,x)+ f N{x,mm (59) 

Jo 

where the kernels to be determined M and N are defined on the triangular domain T. Deriv¬ 


ing (|58]),(|591) with respeet to spaee and time yields the following kernel equations 
for 1 < / < n, 1 < j <m 

n m 

AAMij(x,^)-pjdpMij{x,^) = + J]K~Npj(x ,0 

p=\ 


(60) 


k=\ 


for 1 < / < m, 1 < j <m 


Pid^Nij{x,^) + Pjd^Nij{x,^) = ^o-i^Mkj{x,^) + Npj{x,^) 


(61) 


k=l 


p=l 


April 29, 2015 


DRAFT 














13 


along with the following set of boundary eonditions 


'i\ < i < m, \ < i < n 

Mij{x, x) = 

A... 

(62) 

J: 

-h 

1 

C: 

VI < Uj <m, z y, 

Nij{x, x) = 

0 

(63) 


besides, evaluating ([58]),([591) at jc = 1 yields 


VI < j < i < m 


Ntj(l,x) = Y,PikMkjihx) (64) 

k=i 


To ensure well-posedness of the kernel equations, we add the following artifieial boundary 
eonditions for Nij(i < j) 


VI < i < j < m, Nij{x, 0) = 0 

while the dt, d~j and hjj are given by 

n 

hij(x) = Nij(l,x) - ^pikMtjil,x) 

k=\ 


m px m 

dlj(x, ^) = - ^ Mikix, + I ^ Mikix, s)dlj{s, ^)ds 

k=\ k=l 

m ni 

djj{x, ^) = - ^ Nikix, OcrlJ + I ^ Nik(x, s)d^j(s, ^)ds 

k=\ k=\ 


(65) 


( 66 ) 


(67) 


( 68 ) 


provided the M and N kernels are properly defined. Interestingly, the well-posedness of the 
system of kernel equations of the observer (|60|)-(|65|) is equivalent to that of the eontroller 


kernels s-®. Indeed, eonsidering the following alternate variables 

Mij{x,y) = M,/l -y, 1 -X) = M^j(x,a 

= Nijil -y,l-x) = Nijix,^) 

yields 

for 1 < z < n, 1 < j < m 


(69) 

(70) 


PAMij(x,y) - AidyMij(x,y) = - ^ crl^Mkj(x,y) - ^ Npj(x,y) (71) 

p=i 


k=l 
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for 1 < / < m, 1 < j < m 


f^A^ijOc, y) + f^idyNi j(x,y) = -Y^ Mk j(x, j) - ^ cr,> Npjix, y) (72) 

k=l p=\ 

along with the following set of boundary eonditions 


1 < i < m, 1 < j <n Mij(x,x) = , = mij 

in + dj 

(73) 

VI < z, y < m, i + j, NijOc^X) = 0 

(74) 

n 

VI < y < z < zw, Nijix, 0) = ^ PikMkjix, 0) 

Z:=l 

(75) 

Vl<z<y<zw, Nij{l,y) = 0 

(76) 


whieh has the exaet same strueture as the eontroller kernel system, the well-posedness of whieh 
is assessed in Theorem 13.21 


C. Output feedback controller 

The estimates ean be used in an observer-eontroller seheme to derive an output feedbaek 
law yielding finite-time stability of the zero equilibrium. More preeisely, we have the following 
Lemma. 

Lemma 4.2: Consider the system eomposed of the original Q-Q and target systems (|T7])- 


(|49]) with the following eontrol law 

U(,t)= f [K(l,f)u(^) + L(l,OHf)]d^-^M,l) 
Jo 


(77) 


where K and L are defined by (|3^-(|4T]). Its solutions {u,v,u,v) eonverge in finite time to zero. 


Proof Proposition |4T] along with the existenee of the observer baekstepping transformation 
yields eonvergenee of the observer error states u, v defined by ([50|)-([5^ to zero for t > 


Therefore, for t > tp, one has v{t,0) = v(t,0) and Theorem 3.3 applies to the observer sys 


tern (47H49). Therefore, for t > Itp, one has (u,v,u,v) = 0 whieh also yields {u,v) = 0. 


*the proof of this claim follows the exact same steps as in the controller case, see Section 


III-C 
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V. Motion planning for homodirectional systems 
A. Definition of the motion planning problem 

Consider now the case n = Q. Then system Q-Q reduces to 

vft, X) - h~vfit, X) = S~“v(r, x), 

where coefficients A“ and 1r~ defined as in Q and ([^, with boundary conditions 

v{t, 1) = U{t). 


(78) 


(79) 


For simplicity in this section we drop the super-indices in the coefficients. 


Equation (78) represents a system of m states moving in the same direction (in this case, from 
right to left). We call such a system homodirectional (in oposition with heterodirectional systems, 
whose states move in different directions, such as <13-0 with n,m 4^ 0). Homodirectional 
systems are inherently finite-time stable. Physically, this is due to the fact that they are transport 


equations with information flowing only in one direction; thus, setting IJif) to zero in (79) and 
solving the equations with the method of characteristics, we obtain u{t,x) = 0 for t ^ (the 
slowest transport time in ([78])). 


For (78)-(79) we consider the following motion planning problem. Given 0(0, a known 


function defined as 

O(0 = (oi(0 o„(0J 

find the value of U{t) so that v(r, 0) = 0(0 for t > tM, for some tM > 0. 


(80) 


Remark 7: Even though the plant (78) is finite-time stable, and a formula for the states can 
be written by using the method of characteristics, the motion planning problem is not trivial to 
solve. The entanglement of different states moving with different speeds severely complicates 


finding a solution. This design difficulty will be explicitly shown with an example in Section V-C 


B. Tracking control design 

The following result solves the motion planning problem. 


Theorem 5.1: Consider system (78) with boundary conditions (79), initial condition vo 6 


(L^(0,1))'”, and feedback control law 

i=m 1-1 

^ -V r 'L,, 


um = 'ii\i+ 


HI 






L,,(^,0)O,H + 


1-^ 

di 


d^ (81) 
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Then, v{t,0) = 0(0 if t > tu, for tu = ZJ=i 

Remark 8: The motion planning problem has been solved for the homodireetional case for the 
sake of clarity. However, it can be formulated for the full heterodirectional system ([T])-Q with 
only minor modifications. Noting u{t,0) = Qovit,0), the values of some m,’s could be chosen as 
part of the output instead of some of the v,’s, for a total of m states. The only condition would 
be that all the rows of the output vector (written in terms of the v,’s) are linearly independent. 


Proof: We start by using the backstepping transformation (28)—where the kernels K are 
zero due to n being zero—to map (|78])-(|79|) into the target system 

A(t, - AySx(t, = G{x)m. (82) 


where G{x) was defined in Section III-B as a function of the kernels, with the following boundary 
conditions 


1 ) = 5 ( 0 , 


(83) 


where B{t) in ( [83] ) is a function defined as 

5(0 = (bi( 0 5„(o) , (84) 

with components to be determined. B represents an extra degree of freedom that did not appear 
in the target system for the homodireetional control problem (Equation [T4| ). It will be used to 
solve the motion planning problem. The presence of B{t) in the boundary conditions does not 
change the backstepping transformation; however it modifies the feedback control law to 


/' 


U{t) = B{t)+ 


(85) 

Now, noticing that if one sets x = 0 in the transformation ( [28] ) one obtains v,(t, 0) = Pft, 0), it is 
clear that we only need to solve the motion planning problem for the target (8 system by using 
B{t). The next steps of the proof are devoted to finding the value of B{t). 


Using the method of characteristics, the explicit solution for each state fi{t,x) of (83) with 


boundary condition ([83]) at time t > — is 


Pi{t,x) = Bit+'^\ + - f G(Oj8lt+^,o]d^, 

\ Pi I PiJx \ Pi I 


Using ( [T^ and (42) in (86), we obtain 


/3i{t, x) = Bi\t + 


Pi 


ftX ,r’ 




Pi 


( 86 ) 


(87) 
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To solve now the motion planning problem, consider first (87) for i = 1 and x = 0, for t > ^. 


Imposing /?i(t,0) = Oi(t), we obtain: 




( 88 ) 


thus, setting Bi{t) = Oi for t > 0, we obtain the desired behavior for 0) for ? ^ 


Now consider (87) for i = 2 and x = 0, for t Imposing /32(t,0) = 02(0, we obtain: 


02(0 = ^21? 
Solving for B 2 as before 

B 2 (0 = O 2 1 1 + 


B2I Jo 




-)- f- 

B2I Jo B2 


L2i(^,0)A|t+^,0lJ^. 


(89) 


(90) 


To be able to substitute for Oi(0 in the whole domain of the integral in (90) we need to 

wait until ^ Thus choosing 


B 2 (0 — O 2 1 1 + 


-)- f- 

B2I Jo M2 


L2i(^,0)0ilt+^,01J^, 


(91) 


we get that /32(t,0) = 02(0 for t > ^ ^ (as we have to wait an extra ^ time for (91) to 

propagate). It is clear that this procedure can be continued for / = 3,..., m. Thus we obtain that 


5J0 = O,(t+-)-y r ^L,X^,0)O, (t+^ 


Mi 


^Jo Mi 


Mi 




(92) 


solves the motion problem for /3; for t > li/=i Applying (92) for z = 1,..., m and substituting 
in (|85]) produces the feedback law ([81]), thus solving the motion planning problem in time 


tM 


Z m J_ 
7=1 /i/ 


Remark 9: Theorem 5.1 gives in fact tracking (in finite-time) of the desired output signal, a 
result stronger than pure motion planning. To obtain a pure motion planning result, one should 


take (87)—the explicit solutions of the target system obtained in the proof of the theorem—and 


substitute the values of 5, found in (92), so that the jS,’s are explicit functions of the O/s. Then, 


using the inverse backstepping transformation (45), find the v,’s as explicit functions of the 0,’s 


and substitute them in the control law (81), which would then be an exclusive function of the 


outputs. We omit this result for lack of space. 
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C. An explicit motion planning example 

Next we present an speeifie example of a motion planning problem for m = 2. Consider the 
plant 


with boundary eonditions 


Vit{t,x) - PiVu{t,x) = Cri2V2(t,x), 
V2t(t,x) - P2V2xit,x) = Cr2iVi{t,x), 

V,it, 1) = Uiit), V2(t, 1) = U2. 


(93) 

(94) 

(95) 


The objeetive is to design U\{t) and U 2 {t) so that vi(t,0) = Oi(t) and V 2 (t, 0) = 02(t) for some 


functions 0i,02 for t > tM- Notice that since (93)-(95) is explicitly solvable, one might think 
that the inputs can be directly designed. Using the method of characteristics to explicitly write 
a solution of the system, one gets, after time t = 


vi(t,0) = Ui\t 


V2{t,0) = U2\t 


-h-f 

hi Ml Jo 

-)^-r 

hi Ml Jo 


(TnV2\t- 

Ml 

Cr2iVi |t - ^,^1 


(96) 


(97) 


However, if one tries to proceed as in the proof of Theorem |5.1[ by plugging in Oi(t) in (96) 
and 02(0 in ( |97| ), and then solve for UJt) and U 2 {t), one ends up with a feedback law that 
requires knowing/MtMre values of vi and V 2 , i.e., a non-causal (and therefore not implementable) 
feedback law. Thus, a direct approach does not work even for the m = 2 case. To solve the 


motion planning problem, we resort to Theorem 5.1, in this particular case, the motion planning 
problem is solved by the inputs 


f/i(0 = Oi t+- 

Mi 


ii.r 

h! Jo 


f 


Ln(h^)vjm+ Lnih^)v2{m, 


(98) 


f/2(0 = ®2l^+^)- ^7.2l(^,0)Oi|t+^jj^+ I L2 i(1,^)Vi(^)J^+ I L22{h^)V2(Od^, 

(99) 


f 


f 
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Ln{x,0 


Vo'i20~2i / 2 j cri2(r2iix - 

H2-I^i ^ ^^\{x-^) '\//i-yU2V yUl 

0 , 


£> 

^ m 

^<-x 

^ fii 


(106) 


il2(;c,^) 


trzi j / 2 j cri20-2i(x - - H2^ 

“Ul -^2 V /'I 

0 , 


’ ;ji 


(107) 


i2i(;c,^) 


0-21^^ 

t^lX-H2^ 


f 2 

0-i20-2l(x - 0(MiX - ^2^) 


a-2iMx-0 j 

f 2 

crnO-2i(x - ^){fiix - fi2^) 


^2 

O-I2(MiX-/J20^ 

Ml -M2^ 

M2 


(108) 


i22(;c,^) 


j cri2<r2i j 

2 

10-120-21 (X - ^)(p.iX - IX 2 O 

S M2(x - ^)(MiX - M2^} ' 

Ml -M2 ’’ 

V ^2 J 


(109) 


where the kernels Ln, L 12 , L 21 and L 22 satisfy 

Hid^Ln{x,^) + ^iid^Ln{x,^) = o-2iLi2(;c,^) 
Hid^Lnix^O + l^2d^Lnix,^) = crnLnCx,^), 
H2dxL2l(x,^) + lild^L2l{x,^) = (T2 iL22{x,0 
IJ.2dxL22ix,^) + l^2d^L22{x,^) = (Ti2L2i{x,0, 

with boundary eonditions 


Lii(a:,0) = Li2(;c,0) = L22(x,0) = 0, 


Lx 2 ix, x) = 


cr 12 


L2i(x,x) = 


cr 21 


( 100 ) 

( 101 ) 

( 102 ) 

(103) 

(104) 

(105) 


yU2-jUi yUl-jU2 

plus the artifieial boundary eondition L 2 i(l,^) = / 2 i(^), where the funetion hi is arbitrary. These 
kernel PDEs ean be explieitly solved using teehniques akin to those used in [[25l . The resulting 


kernels (whose validity ean be verified by substitution in the kernel equations) are given by ( |106[ )- 
( |109| ). where lo and Ii are the modified Bessel funetions of order 0 and 1, and Jo and Ji are the 
(regular) Bessel funetions of order 0 and 1, respeetively. 

The kernels appearing in (|9^-(|97|) are depieted in Fig for the ease = 1, /12 = 0.2 and 
cri 2 = 2, 0-21 = 5. It ean be seen that Ln(l,^) and Li 2 (l,^) have a monotone behaviour (they 
are always negative or zero), whereas L 2 i(l,^), L 2 i(^, 0), and L 22 (l,^) are oseillatory. Fig. 
shows Lii and L 12 in the whole domain 7“; notiee that Li 2 {x,^) is diseontinuous along the line 
^ “ /o ('''l^l‘-h is the lower domain on Figure whereas Lii(ji;,^) is not diseontinuous. On the 
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Other hand, it is evident that hiiO = L 2 i(l,^) is rather non-trivial. In fact, the procedure that 
was followed to find these explicit solutions was not setting a value of hi a priori, but rather 
extending the domain shown in Figure ^up to x = so that boundary condition (105) can 
be used to actually find the value of hi- 



Fig. 1: Motion planning kernels (n = 0, m = 2). Solid: Lii(l,^) and L 2 i(^, 0). Dash-dotted: 
Li 2(1,^). Dotted: L 2 i(l,^). Dashed: L 22 (l,^). 


LiAx,^) Li2(x,^) 



Fig. 2: Motion planning kernels Lii(x,^) and Lnix,^) (n = 0, m = 2). 
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VI. Proof of Theorem [3]^ well-posedness of the kernel equations 

To prove well-posedness of the kernel equations, we elassieally transform them into integral 
equations and use the method of sueeessive approximations. 

Remark 10: Similar proofs have been derived for less general systems, e.g. in [|3 or ITT^ . 
The proof is more involved here due to the existenee of homodirectional controlled states, whieh 


lead to the homodireetional kernel PDEs (37). 


A. Method of characteristics 

1) Characteristics of the K kernels: For eaeh l<i<m, l<j<n, and (x, 6 T, we define 


the following eharaeteristie lines eorresponding to Equations ( [3^ 

-(x, s) = -fii, s 6 [O, sffx, f)] 


‘ dxij 


ds 


Xij(x,f; 0 ) = X, Xijix,^; sUx,^)) = xUx,^) 


( 110 ) 


ds 


-(x,^; s) = Aj, s e [o, 


( 111 ) 


^,/(x,^;0) = fiAx,^;sfAx,f)) = xfAx,f) 


ds 


These lines, depieted on Figure originate at the point (x, f) and terminate on the hypothenuse 
at the point (x^.(x,^),x^(x,^)j. The expressions of Xij{x,f; 5 ), fifx,^; s) sfj(x,f) and x^(x,^) are 
omitted here beeause of laek of spaee, but are straightforward to eompute. Integrating ( [361 ) along 

these eharaeteristie lines and plugging in the boundary eondition ( [38] ) yields 

X sfjix^) r n m 

^ crljKikiXijix,^; s),fij(x,^; s)) + ^ o-~p]Lip{xij{x,^-, s),fij(x,^; s)) 

-/i=i p=i 

( 112 ) 

2 ) Characteristics of the L kernels: For eaeh 1 < / < m, 1 < 7 < m, and (x, f) 6 T, we define 
the following eharaeteristie lines {Xij{x,f\-),t^ij{x,^\-)) eorresponding to Equations ( [37] ) 

(x, v) = djiii, V 6 [ 0 , vj(x, ^)1 

(113) 


dv 


Xij(x,^; 0 ) = X, XiM,f;vfj(x,f)) =^J.(x,^) 

^(x, y) = Cijiaj, v 6 [o, vj.(x, f)\ 
^ij(x,f-0) = f, ^ij(x,f;vfAx,f)) = d(x,0 


(114) 
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Fig. 3: Characteristic lines of the K kernels 


where 6,, is defined by 




1 if z > j 
-1 otherwise 


(115) 


These lines all originate at (x,^) and terminate on ffT at the point They 

are depicted on Figures |4]^ in the three distinct cases i < j, i = j and i > j. The detailed 
expressions of;^f,•y(x,^; s), ^ij(x,^; s) vfj{x,^),xfj{x,^) and ^[j{x,^) are, again, omitted here because 
of space constraints. Integrating ( |37] ) along these characteristics and plugging in the boundary 
conditions (|^,(|^ and yields 


1 ” 

Lij (x, = 6 ij(x, ^)lij + (l - Sij{x, ^)) — ^ ArqrjKirOcfjix, ^), 0) 




X vfjixS m n 

_p=i k=i 


dv (116) 
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Fig. 4: Characteristic lines of the kernels L,y for i > j 


where the coefficient 6ij(x, ^), defined by 

0 if / = j 

= < 0 if / < 7 and - ^jX < 0 ^ (117) 

1 otherwise 

reflects the fact that some characteristics terminate on the ^ = 0 boundary of T, while others 


terminate on the hypotenuse or on the x = \ boundary of T. Plugging in (112) evaluated 
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Fig. 5: Characteristic lines of the kernels L„ 


at (Xu{x,^),0) yields 


1 ” 

Lij (x, = 6ij{x, ^)lij + (l - Sij(x, ^)) — ^ Arqrjh 




X 1 ■A " 

+ (1 - Sij(x, ^)) — y Arqrj I y o-l;!'Kik(xir(xfj(x, ^), 0; s), ^iriXijix, 0, 0; s)) 

r=l "^0 Li=l 


^ a-p^ Lip{xir(xfj{x, ^), 0; 5 ), ^trOcfjix, ^), 0; 5 )) 


p=i 


J 5 


X vfjixS m n 

Yu ^~P~i ^ Zi 

_p=l j(:=l 


Jv (118) 
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Fig. 6: Characteristic lines of the kernels L,y for i < j 


B. Method of successive approximations 


We now use the method of successive approximations to solve equations (112),(118). Define 


first 


VI < z < m, 1 < 7 < R, (pij{x,f) = kij, (119) 

VI < z < zn, 1 < j < m, 

1 ” 

fijix, = 6ij(x, ^)lij + (1 - Sijix, ^)) — X '^'■qrjkir (120) 

r=l 
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Besides, we define H as the vector containing all the kernels, reordered line by line and stacked 


up, and similarly (p, as follows 


Hi 


'Ku 


( \ 
01 


( \ 
<Pn 

Hfim 


Kfnn 

, 4> = 

4^nm 


^mn 

Hnm+i 


Lu 


4^nm+l 


*Aii 

H 7 




^4^nm+nP- ^ 




( 121 ) 


We consider the following linear operators acting on H, for 1 < / < m, 1 < j < n 
0 ,y[H](v:, ^) = I ^ cr^lKik(xij(x, 5 ), ^ij(x, s)) + ^ crj Lip(xij(x, 5 ), ^ij(x, s)) 

^0 Vk=\ 0=1 


ds 
( 122 ) 


and for 1 < / < m, 1 < j < m 
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Define then the following sequence 


HVa = 0, (124) 

= <p{x,0 + (125) 


^n(x,^) + On[H'?-'](x,^) 
<^i„(^,^) + Oi„[H'?-'](z,^) 

^2i(^,^) + 02i[H^-'](^,^) 

ifinnix,^) + 0„,„[H‘?“^](;c,^) 

ilfn(x,0 + ^n[^^-'](x,0 


(126) 


One should notice that if the limit exists, then H = lim H^(x, ^) is a solution of the integral 

q—^+oo 

equations, and thus solves the original hyperbolic system. Besides, define for q > I the increment 
AH^ = with AH° = ^ by definition. Since the functional O is linear, the following 

equation AH^(x, ^) = 0[H^"'](x, ^) holds. Using the definition of AH^, it follows that if the sum 

+ 00 

2 AH^(x,^) is finite, then 

q=0 


+ 00 

H(x,^) = ^AH^(x,^) (127) 

q=0 

In the next section, we prove convergence of the series in L“. 


C. Convergence of the successive approximation series 

To prove convergence of the series, we look for a recursive upper bound, similarly to, e.g. [|l2l- 
More precisely, let e be such that 


0 < e < 1 - max 


(ih 


l<7</<m 


(128) 


Then, the following result holds 

Proposition 6.1: For q > I, assume that 

V(x,^) 6 7“, Vi = 1,..., nm + 


\AHfx,f)\ < 


M‘i(x-{\-e)^y 
f -i- 

q\ 


(129) 
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then, it follows that 

V(x,^) 6 T, V/ = l,...,m, vy = l,...,n, |0,y[AH](x,^)| < 0- (130) 


and 


V(a:,^) 6 T, Vi = l,...,m, Vj = l,...,m, 


(V+1)! 


( 131 ) 


(«+!)! 

The proof of this proposition relies on the following Lemma, whieh is erueial and different with 
previous works. 

Lemma 6.2: For q £ N, (x,^) £ T, and sfj(x,^), vfj(x,^), x,-j(x,^,-), ^i/x,^,-), Jtfqix,^,-), 
^ij(x,^,-) defined as in ( |110| ),( [T]T| ),( |113| ),( ri 14| ), respeetively, the following inequalities holds 
VI < i < m, VI < j < n 

J [Xij(x, ^;s)-(l- e)^ij(x, ^;s)) ds < M^- 


VI < i, j < m 


X yfjLS 

v) - (1 - eXij(x, v)f dv < M, 


where 


Mi = max 


1 


q + \ 

{x-{\-e)0^^^ 
q + \ 

1 


(132) 


+ (1 - e)Aj _ (1 _ 

Proof: Consider the following ehange of variables, noting ( |110| ),( [TTT] ), 

T = Xij(x,f; 5) - (1 - 6)fij{x,f; s). 


dr 


dxii d^ii 

5)-(l -6)^(x,f; s) 
ds ds 


ds 


{-Hi - (1 - e)A^ds 


The left-hand-side of ( |132| ) beeomes 

px".\x. 

{xij(x, 5) - (1 - 6)fij{x, f-, s)Y ds = 

0 Jx-(1-. 








Hi + (1 - e)A 


■dr 


(133) 

(134) 

(135) 

(136) 

(137) 

(138) 


(X - (1 - e)fr^ - (xfjix,f) - (1 - e)f[j(x,0) 


(139) 

?+i 


< M. 


{Hi + {l-e)Aj){q+l) 

q + \ 


(140) 

(141) 
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where we have used the faet that for all 1 < / < m, 1 < j < n, one has 

( 142 ) 

whieh is trivially satisfied sinee (xfj(x,^),^[j(x,^)) e dT and e > 0. Consider now the following 
ehange of variables 


T 5) - (1 - e)^ijix,^; s), 

dT 


dXa d^a 

5) - (1 - 6)^{x,^; s) 
ds ds 


ds 


= Qj (/i; - (1 - e)/i;) ds 


(143) 

(144) 

(145) 


Thus, the left-hand-side of ( |133| ) beeomes 

, y) - (1 - eXijix, v)f dv = f 

Jx- 






-dT 


x-(i-6)f eij(^iUi - {I - e)iUj^ 

ix-{l- e)^r^ - (xfjix,^) - (1 - eXfjix,^) 


(146) 

<?+! 


; (//,• - (1 - 6)/2y) (q 1) 


(147) 




(148) 


Given the definition of 6,y given by ( |115[ ), one has 

|yU, -(1 -6)yUy if / < 7 
[{1 - e)jTj - lUi if i>j 

Therefore, given the definition of e (Equation ( |128| )) in the ease i > j and the ordering of the //, 
in the ease i < j, one has 

-6ij (fXi - (1 - e)yUy) >0 (149) 

Besides, sinee (Xn(x,^),^f.(x,^)) 6 T, one has “ (1 “ > 0 and (|147|) beeomes 

u V V V U / ' ' 


J (xij(x, v) - (1 - e)^ij(x, y)) dv < Ma 


{x-{l-e)^r^ 

q -i- \ 


(150) 


whieh eoneludes the proof. ■ 

Remark 11: Notiee that ( |149[ ) also implies that, for any (x, ^) 6 T and 1 < / < m, 1 < 7 < n 
the funetion 


y 6 [0, v’-j{x, ^)] Xijix y) - (1 - e)^,7(-^’ 


(151) 
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is strictly decreasing, in particular the following inequality holds 


0 <- (1 - <x-{l-e)^ 


( 152 ) 


whieh will be useful in the proof of Proposition |6.1[ 


p-|. 

Ai J 


Proof of Proposition 6.1 ■ Define 

d = max{/l„,/ii}, d = max- 

a = max (t~^, a~~], q = max{^,y} 

>J ij 

M = (nAAq + l) (n + m)&Mx, 

0 = max max |0,j(x,^)|} 


Let now 6 N and assume that 

^{x,^) 6 T, Vz = 1, ...,nm + 

Then, for 1 < z < m, 1 < j < n, {x, f) &T one has 

pfjixS n 

IO; / [AH] (x, I < I y aljAKikiXij(x,^;s), ftfx, 

^0 Z:=l 


\AHi(x,0\ < 0 -;- 


(153) 

(154) 

(155) 

(156) 

(157) 


^)) 


+ Zz ALip(xij{x, s), fij{x, s)) 

p=\ 


ds (158) 


using ( |132| ) and ( |157| ), this yields 

|0,y[AH](x,^)| < in + m)d-- f 

Jo 


fM) _M^(x,j(x,f; s)-a- e)f^j(x,e, s)f 
^ -A- 


ds 


. . , (^-(i-^)r' 

< (zz + m)cr - M^- 


^ + 1 


< 0 


M‘i^\x - (1 - 6 )^)?+' 

(^+ 1 )! 


(159) 

(160) 
(161) 


Similarly, for 1 < z, j < m, one gets, using ( |157[ ) 

I ¥,7 [AH] (x,^) I < AAqin+m)d-Y\ f 

TTT-'o 

+ (n + /n)cr 


(xp,ao) _M^(x,,(xfj(x,aO;s)-(l-e)f^,(xfj(x,aO;s)^^ 

0 - 


ql 


J 

Jo 


tef) _ (xifx, eXijix, f; v)f 




q\ 


dv (162) 


ds 
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Then, using ( |132| ) at (x,^) = (xfi(x,^),0) and ( |133| ) yields 


|'P,y[AH](x,^)| < AAq(n + m)d-n^M^M‘< 


(-?+!)! 


, ^_-MmAx-{l-e)^r^ 

+ (n + m)(T(p -;-—- (163) 


Inequality ( |152| ) yields 

|'F;j[AH](;c,^)| < [nAAq + l) (n + m)d-^M^ 


{q+l)\ 

M\x - (1 - 

(^+1)! 


< 0 - 


(^+ 1 )! 


(164) 

(165) 


which concludes the proof. 


Proposition 6.1 directly leads to Theorem 3.2, since by the same procedures presented in DTI 
and [fT^ . one has that ( |127[ ) converges and 


|H(ji;,^)| 


^AH^(x,^) 

q=Q 


< (pe 




(166) 


VII. Concluding remarks 

We have presented boundary control designs for a general class of linear first-order hyperbolic 
systems: an output-feedback law for stabilization of heterodirectional systems and a tracking 
controller for motion planning for homodirectional systems. 

These results bridge the gap with the results of, e.g. where the null (or weak) control¬ 
lability of (n -I- m)-state heterodirectional states is proved but no explicit design is given. 

Our results open the door for a large number of related problems to be solved, e.g. collocated 
observer design, disturbance rejection, similarly to [[H, parameter identification as in [fTTlI . output- 
feedback adaptive control as in Q, and stabilization of quasilinear systems as in [ITII. 

Another important question concerns the degree of freedom given by Equation ( |4T] ) in the 
control design. The effect of the boundary value of the kernels on the transient performances of 
the closed-loop system is non-trivial, yet crucical for applications. 
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